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Coupling of Josephson-phase and spin-waves is theoretically studied in a superconduc- 
tor/ferromagnetic insulator/superconductor (S/FI/S) junction. Electromagnetic (EM) field 
inside the junction and the Josephson current coupled with spin-waves in FI are calculated 
by combining Maxwell and Landau-Lifshitz-Gilbert equations. In the S/FI/S junction, it 
is found that the current-voltage (I-V) characteristic shows two resonant peaks. Voltages 
at the resonant peaks are obtained as a function of the normal modes of EM field, which 
indicates a composite excitation of the EM field and spin-waves in the S/FI/S junction. We 
also examine another type of junction, in which a nonmagnetic insulator (I) is located at one 
of interfaces between S and FI. In such a S/I/FI/S junction, three resonant peaks appear in 
the I-V curve, since the Josephson-phase couples to the EM field in the I layer. 

KEYWORDS: Fiske resonance, Josephson junction, superconductor, ferromagnetic insulator, 
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1. Introduction 

The dc Josephson effect is characterized by the zero-voltage current through a thin in- 
sulating barrier sandwiched by two superconductors. 1 ) This effect is a macroscopic quan- 
tum phenomenon involving phase coherence between two superconductors. When a finite 
voltage(V)-drop appears in the junction, the difference in the phase of superconducting order 
parameter, i.e. Josephson-phase (9), oscillates with time according to 86/ dt = (2e/h)V, and 
the alternating current with frequency [2e/H)V flows in the junction. This ac Josephson effect 
is derived by the gauge invariance including 9. The electromagnetic response dominated by 9 
shows a resonant behavior in the junction. When a dc magnetic field and the dc voltage are 
applied to the junction, the electromagnetic (EM) field is generated by spatially modulated 
ac Josephson current. In this case, the current-voltage [I-V) curve exhibits resonant peaks 
due to the resonance between the ac Josephson current and the EM field generated by the 
spatially modulated ac Josephson current itself. This is called Fiske resonance. 2 ^ 

In recent years, a ferromagnetic Josephson junction composed of ferromagnetic metal (F) 
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and superconductors (S's), i.e., S/F/S junction, has received much attention. 7 ~ 10 ) One of the 
interesting effects is the formation of tt state arising from the Zeeman splitting in F. In ad- 
dition, the interaction between Cooper pairs and spin waves in F is also of importance in 
the transport properties in the S/F/S junction. n ~ 19 ) In a small junction, where the junction 
width is smaller than the Josephson penetration depth, the spin-wave excitation induced by 
the ac Josephson effect is observed. 18 - ) In the recent experiment in a S/F/S junction includ- 
ing a nonmagnetic insulator (I) in one of interfaces between S and F, it has been reported 
that the Fiske resonance has multiple structures that must be associated with the spin-wave 
excitation. 20 ) Volkov et al. have theoretically studied collective excitations in such a junc- 
tion and reported an additional structure in the Fiske resonance induced by spin- waves. 19 ) In 
their theory, a nonmagnetic insulator is crucial to obtain the Fiske resonance coupled with 
spin-waves. On the other hand, another type of ferromagnetic Josephson junction composed 
of ferromagnetic insulator (FI) and two S's, i.e., S/FI/S junction, is also expected to show the 
similar multiple structures in the Fiske resonance. It has been reported that the dissipation 
effect in the S/FI/S junction is smaller than that in the S/F/S junction. 21, 22 ** Such a small 
dissipation in the S/FI/S junction is due to the small probability of quasi-particle excitation in 
the FI. 21,22 ) The damping of spin- waves induced by the similar mechanism is also very small 
in the FI compared to the case in F. 23,24 ) Therefore, the coupling between Josephson-phase 
and spin-waves can be observed more clearly in the S/FI/S junction. 

In this paper, we theoretically study a composite excitation of the Josephson-phase 
and spin-waves in the S/FI/S and S/I/FI/S junctions. First, we calculate the dynamics 
of Josephson-phase coupled with spin waves by using Maxwell and Landau-Lifshitz-Gilbert 
(LLG) equations. Second, we derive the dc Josephson current induced by the Fiske resonance. 
In the S/FI/S junction, two resonant peaks appear in a current-voltage curve for each mode 
of the EM field. These two resonant peaks may be associated with the direct coupling between 
spin-waves and the EM field inside the junction. We also discuss the Fiske resonance in the 
S/I/FI/S junction. The non-magnetic high resistive layer is sometimes important, since the 
magnetic dead layer exist in the ferromagnetic insulator near the S/FI interface. Our results 
clearly show the difference between S/FI/S and S/I/FI/S junctions in the dispersion relations 
of the Fiske resonance. In such a S/I/FI/S junction, we show that three resonant peaks appear 
in the I-V curve for each mode of the EM field. 

The rest of this paper is organized as follows. In Sec. II, by combining the Maxwell and 
LLG equations in a S/FI/S junction, we formulate the dc Josephson current induced by the 
Fiske resonance. In Sec. Ill, the Fiske resonance is discussed in S/FI/S and S/I/FI/S junctions. 
Summary is given in Sec. IV. 
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2. Formulation of Fiske resonance in S/FI/S junction 

The system considered is a Josephson junction with a FI sandwiched by two s-wave 
superconductors (S's) as shown in Fig. 1. The magnetization in the FI is parallel to the 
z-direction. 18 ' A uniform dc magnetic field is applied in the x-direction. In the measurement 
of the Fiske resonance, the dc magnetic field is smaller than several tens of gauss. Therefore, 
we can neglect the in-plane magnetization induced by the applied dc magnetic field. Here, we 
consider that the ac electric and magnetic fields are in the z- and x-direction respectively, both 
of which are uniform in the x-direction. We consider the situation, in which the z-dependence 
of the electric and magnetic fields in the FI is negligible due to the very thin thickness of 
the FI (dpi). In the S regions, it is assumed that the magnetic field depends on y- and z- 
component. The current density has a nonzero y-component in the superconducting regions 
(Meissner current) and a nonzero z-component in the ferromagnetic region (quasi-particle and 
Josephson currents). Based on the above assumptions, the Maxwell equation in each region 
is given by 

d 

rot[E z (y,t)e z ] = - — [n H x (y, z,t) + M x (y,t)]e x , (1) 
rot[H x (y,z,t)e x ] = Jl l (y,t)e y , (2) 

d_ 

dt 



d 

rot [H x (y, t)e x ] - — [D z (y, t)e z ] = J^{y, t)e z + J^{y, t)e z , (3) 



M x (y, t) = / dy'dt' X x(y ~y',t- t')H x (y', t'), (4) 



JI(y,t) = J cS m9(y,t), (5) 

J^y,t) = j^E z (y,t). (6) 

Here, (i = x,y,z) is a unit vector, E z (y,t) is the electric field in the FI, H x (y,z,t) 
and H x (y,t) are the magnetic fields in the S and the FI, respectively. The electrical flux 
density, D z (y,t), in FI is given by D z (y,t) = epiE z (y,t), where epi is the dielectric constant 
in FI. Jj(y,t) and J c are the Josephson current and the Josephson critical current densities, 
respectively. Jq(y, t) and Rp\ are the quasi-particle current density and the resistivity of FI, 
respectively. The x-component of the magnetization, M x (y,t), in the FI is given by Eq. (4). 
The motion of magnetization is described by the Landau-Lifshitz-Gilbert (LLG) equation, 23 ) 

dM 



dM TT a 

= - 7 M x H cS + — 



M x 



dt 



(7) 



dt 1 M 
where M is the magnetization of FI, 7 is the gyromagnetic ratio, and a is the Gilbert damping. 
The effective field, to which M responds, is given by H c g. 

By using Maxwell and LLG equations, we can obtain the voltage coupled with spin- waves 
(the detail of derivation for Eq. (8) is given in Appendix A.) as follows: 



d 2 V(y,t) 



d 2 V(y,t) d FI 1 r° , . , ,,d 2 V(y',t>) 

2 — ^2 h 1 — / dy dtxiy -y ,t-t^ 
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Fig. 1. (Color online) Schematic figure of a S/FI/S junction with ferromagnetic insulator (FI) between 
two superconductors (S's). dyi is the thickness of FI. L and W are the widths of the junction. M 
and H ex are the magnetization in the FI and the applied dc magnetic field, respectively. 



dV(y,t) dpi 1 r° ,,dV(y',t>) 
+ T F i — 7TT ^ r Fii -ttt / dydtx{y-y ,t-t)- 



oo 
oo 



+ TTT J J^' t ) + x2 T ox, 77 / dy'dt'x{y- y ',t-t')JI{y',t'), (8) 



<9t dpi + 2Al Mq y_oo ' dt 

1 1 dpi 1 

\Jj c Jj{ *' t)+ \p c d F1 + 2\ h ]T 

where c F i = \AWPfi + 2Aj)eFi/U ], Aj = v^/Pe^o^Fi + 2Aj) J c ], and T F i = (e F i-RFi) _1 
are the effective velocity of light in the FI, the Josephson penetration depth, and the damping 
factor caused by quasi-particle resistivity, -Rfi, i n the FI, respectively. 
We look for the solution of Eq. (8) in the form 

V(y,t) = V + v(y,t), (9) 

where Vq and v(y,t) are the dc bias voltage and ac voltage induced by the ac Josephson 
current, respectively. In this case, the phase difference, 6(y,t), between two S's is given by 

9{y,t)=Lo J t-k H y + e 1 {y,t), (10) 

where wj = (2e/h)Vo is the Josephson frequency, &h = 2-7r / uo<iFi-£^ex/3 ) o depends on the exter- 
nal magnetic field, H cx , and <3?o is the magnetic flux quantum. 9\{y,t) is related to v(y,t) by 
the equation, 

Substituting Eq. (9) and Eq. (11) into Eq. (8), we obtain the equation for 9i(y,t) as follows: 



&Qi{y,t) = i 

dy 2 c|j 



d 2 0i (y,t) d Fl 1 [°° , , , . POi (y>,t>) 
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dt d F i + 2A L 7-oo 

1 1 rl i r 00 

+ W^' t) + ?7 rf +9A 7T / <h/Mx{v-i/,t-i!)J!(i/J). (12) 

AjJ c Aj J c CtFI + ^Al ^0 ./-oo 

We expand 9\(y,t) in terms of the normal modes of the electromagnetic field generated by 
the ac Josephson current, 



0!(y,t) = Im 



(13) 



J^e^'cos^y) 

.n=l 

where is a complex number and k n = nir/L. This equation of 9i(y, t) satisfies [d9\/dy\ y =Q = 
[d9\/dy\ y= L = 0, which corresponds to the open-ended boundary condition for v(y,t). We 
consider 9\(y,t) to be a small perturbation and solve Eq. (13) by taking Jj(y,t) to be 
J c sin(a;jt — kny). Substituting Eq. (13) into Eq. (12), g n becomes (see Appendix B ) 

C FI n \ B n — iC n 

9n = —T- fl>{kH,-U>j)— 2 — r, (14 

Aj lo^ - n(k n , -wjjwj + zT F nj(k n , -wj) 



2 f L 

Bn= L J dy C ° S ^ nV ^ C ° S ' ^ 

2 /" L 

Cn= L J dy C ° S ^ n2/ ) Sin ' ^ 
-^j) = 1 + Xx(9, -wj)rfFi/[(dFi + 2Aj)// ], (17) 

where oj n = (cfitt / L)n, and q means kn or k n . In the linearized LLG equation, the magnetic 
susceptibility in the FI is given by (see Appendix C) 

Xx{q,uj)=-jM z —2 — — — 2 2 • (18 

iZg — (1 + or jWj + izattscoj 

Here, Sis is spin wave frequency whose dispersion relation is given by 

tts = ^B + |g 2 , (19) 

where fie = j(Hk~ M z /ho). Hk and 77 are the anisotropic field and the stiffness of spin waves 
in the FI, respectively. 

Next, we calculate the dc Josephson current coupled with spin waves as a function of 
the dc voltage and of the external magnetic field. The function, sin(a;jt — kny + 9i(y,t)), is 
expanded with respect to 9±(y,t) and the dc Josephson current is given by 

1 f T 1 f L 

J dc « lim — / dt- / dyJ c cos(ujt - k u y)9i(y,t). (20) 
T^oo 1 J L J 

Introducing Eqs. (13) and (14) into Eq. (20), the analytic formula of the dc Josephson current 
is obtained as, 

t 2 oo 

J *c = ^E*^)' ( 21 ) 

J n=l 

* n = Re[M(fcH,wj)X], (22) 
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Fig. 2. (Color online) Dc Josephson current density, Jd c , as a function of dc voltage, Vq, in a S/FI/S 
junction. The solid line is the total dc Josephson current. Red, blue, green, purple, and light blue 
lines are the dc Josephson current of each mode number, n, of electromagnetic held. The applied 
dc magnetic field determines n via Eq. (25). 



X — (23) 

w 2 - fi'(k n ,uj)ujj + fi"(k n ,uj)T F iLJj + i [fi'(k n ,ujj)T F iu;j + n"(k n , uj^uj] ' 

A*(9, wj) = fjf(q, wj) + ifi"(q, wj), (24) 
2(f) sin (n<j) — nir/2) 
> + n/2 TT(f) — n7r/2 

where eft is equal to $/$o and $ = fioH^dpiL. fi'(q, coj) = Re[fj,(q,coj)] and fjf'(q,u)j) 
Im[n(q, wj)] (See Eqs. (17) and (18)). 



K 2 ( 



(25) 



3. Results and discussion 

In this section, we examine the numerical solution for Eq. (21). Figure 2 shows the dc 
Josephson current density induced by the Fiske resonance as a function of the dc voltage 
for Aj/L = 1, ^b/^l = 3, v/(^l) = 3 x l(T 16 m 2 , 7 M 2 /(/x wl) = 1, 25) a = 1 X 1(T 4 , 24 ) 
dpi /(dpi + 2Al) = 0.1, Tfi/wl = 3 x 1CP 1 , 26 ) and wl = cfivt/L. is fixed as re/2 in the 
F n ((j)) function. In the solid line of Fig. 2, the normalized dc Josephson current density J^ c 
in Eq. (21) is shown as a function of normalized dc voltage. 27 - > Red, blue, green, purple, and 
light blue lines are the dc Josephson current of each mode number of electromagnetic field, 
re. In this Figure, the resonant behavior of J& c is due to the Fiske resonance in the S/FI/S 
junction. However, at Vo/(oJi l h/2e) ~ 3, it is found that additional structures of Jd c appear. 
From Fig. 2, it is found that two resonant peaks appear for each n. The appearance of the 
two resonant peaks in the S/FI/S junction are very different from conventional Josephson 
junctions, in which a single resonant peak appears for the each mode of EM field. 2-6 ^ Large 
peak around Vo/(uJi,h/2e) = 3 is due to the summation of n > 8 because dc Josephson currents 
of contribution from large n appear around Vo/(wL^/2e) = 3 in a manner to be described. 

To elucidate the origin of the resonant structures in the S/FI/S junction, we analyze 
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Fig. 3. (Color online) Dc voltage, Vq, as a function of mode number, n, of electromagnetic field in 
S/FI/S junction. 



Eqs. (21) and (22). When the denominator of ty n in Eq. (22) is minimum with respect to 
ojj, takes a maximum, so that the dc Josephson current shows the resonant behavior as 
shown in Fig. 2. The dc voltage at which the resonance occurs is determined by neglecting 
the damping term in Eq. (22) as a = Tpi = 0. Setting the denominator of \P n to be zero, the 
voltage is given by 



V 



± 



~2e 



2 , r>2 , d FI jM z Q S , if 2 , n2 , d FI T^^S V , 2 r>2 

w£ + + ; ±\ [u)£ + ni + : iui Sis 

s dpi + 2A L mo V V d F i + 2A L /i / n s 



(26) 



where w n is the frequency of the EM field in the FI, and Og and "fM z {ls/ no are the frequency 
of spin-waves and the real part of the magnetic susceptibility with a = in the FI. We have 
two dc voltages, and Vq - , at which the Fiske resonance occurs for each n. From the analytic 
formula in Eq. (26), it is found that two dispersions result from the coupling between the EM 
field and spin waves in the FI. Figure 3 shows a Vq-u curves obtained by Eq. (26). The vertical 
axis is the dc voltage normalized by uj-^h^e and the horizontal axis is the mode number of 
EM field. In Fig. 3, Vq + and Vq are shown by open circles and open squares, respectively. 
For n < 3, Vq + is nearly constant as a function of n, whereas Vq - is linear with n. In Vq + for 
n < 3, the voltage is nearly equal to Og/i/2e, which relates to the spin wave energy in Fig. 3. 
For n > 3, V + increases as a function of n, whereas Vq~ becomes flat with increasing n. For 
n > 3, V " is nearly equal to f2s^/2e in Fig. 3. Therefore, it is found that the flat behavior 
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Fig. 4. (Color online) (a) Dc Josephson current density, Jd c , as a function of dc voltage, Vq, by chang- 
ing Tfi/wl in a S/FI/S junction, (b) Jdc as a function of Vb, by changing a in a S/FI/S junction. 
Where parameters used by numerical calculation are Aj/L = 1, fie/^L = 3, "fM z /(fj,QUJi J ) = 1, 
d F1 /(d F1 + 2A L ) = 0.1. 



of the voltage comes from the spin- wave excitation in the FI. The spin-wave excitation in the 
FI is induced by the EM field generated by the ac Josephson current and the effect of the 
spin-wave excitation is reflected in the Fiske resonance in the S/FI/S junction. 

Here, we examine the Tpi- and ev-dependence of Jd c . Figure 4 (a) shows J^ c as a function 
of Vq by changing Tfi/wl- From this figure, it is found that the Fiske resonance without 
very sharp peaks around Vo/(uj^h/2e) = 3 exhibits strong damping by increasing Tfi/wl- 
On the other hand, very sharp peaks around V^j {uj\ J fij2e) = 3 almost does not depends on 
Tfi/wl unlike another resonant peaks. Next, we focus on sharp peaks around Vo/(u}i,h/2e) = 3. 
Figure 4 (b) shows Jh c as a function of Vq by changing a. These peaks around Vq/ (u\ J h/2e) = 3 
in Fig. 4 (b) strongly depend on a because these resonant peaks mainly comes from spin-waves. 
From Fig. 4, we can easily obtain the Fiske resonance coupled with spin-waves in the S/FI/S 
junction due to the small T-pi and a. 

The effect of spin-waves having a finite wave number q is neglected in the Fiske resonance 
because of the following reason: In Eq. (19), the first term fie is caused by the anisotropic 
and demagnetizing fields and finite wave number q is given by nn/L. In a conventional FI, 
M7b is about tens of //eV. 23 ) On the other hand, rjq 2 is of the order peV due to the small 
stiffness of spin-waves 28 ^ when the width (L) of the junction is a few mm. 

Next, we consider the Fiske resonance in the S/I/FI/S junction. The details of the calcu- 
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Fig. 5. (Color online) Dc voltage, Vo, as a function of mode number, n, of electromagnetic field 
in S/I/FI/S junction. We take parameters as fie/wL = 3, c\/cpi = 0.5, 'yM z /(/j,oU>i,) = 1 and 
dml(d m + A L ) = 0.2. 



lation are given by Appendix D. To analyze the origin of the Fiske resonance accompanied 
by the spin-wave excitation in the S/I/FI/S junction, we examine Eqs. (D-15) and (D-17). 
The condition of the resonance is given by the minimum in the denominator. To find out the 
voltage at which the Fiske resonance occurs, we neglect the damping term in the denominator 
of Eqs. (D-15) and (D-17). As a result, the voltages at the Fiske resonance are given by 



2e 



h 

2e 



-0J r 



(27) 



d F i jMzth 



«fi + Al ho 



(28) 



indicating that there are three dc voltages resulting in the Fiske resonance for each n in the 
S/I/FI/S junction. Figure 5 shows the Vq-u characteristic obtained by Eqs. (27) and (28). 
The vertical axis is the dc voltage normalized by co^h/2e and the horizontal axis is the mode 
number of EM field. In Fig. 5, V , V + and are shown by open triangles, open circles 
and open squares, respectively. V^ 1 linearly increases as a function of n as shown in Fig. 5 
and comes from the resonance between the EM field in the I and the ac Josephson current. 
For n < 3, V + is nearly flat as a function of n, whereas V ~~ shows linear behavior with n. 
For n > 3, V + increases as a function of n, whereas V ~ becomes flat with increasing n. The 
voltage of flat region as a function of n is nearly equal to ^s^/2e in Fig. 5. 

In this paper , we assumed that the magnetization of FI is a single domain structure. 
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Usually, magnetization structure in the FI possesses a complicated domain structure. The 
ferromagnet with multidomain structure has multiple magnetic resonance modes different from 
the single domain system. 29 ) Therefore, we can expect several additional Fiske resonance peaks 
in the I-V curve arising from the domain structure in the FI. Those domain structure and their 
spin dynamics, which depend on our choice of material, have rich variety and wide spectrum. 
However, it is difficult to include the multidomain structure of FI in the present theory. The 
Fiske resonance in the ferromagnetic Josephson junction with the magnetic domains is beyond 
the scope of the present paper, and will be studied in another one. 

4. Summary 

We have theoretically studied the coupling of Josephson-phase and spin-waves in the 
S/FI/S junction. The dc Josephson current induced by the Fiske resonance is calculated by 
combining the Maxwell and the Landau-Lifshitz-Gilbert (LLG) equations. The dc Josephson 
current shows the resonant behavior as a function of the applied dc voltage. We derived the 
analytic formula of the resonant condition in the Fiske resonance in the S/FI/S junction. Two 
resonant peaks appear in the I-V curve for each mode number of EM field. We found that 
the two resonant peaks are generated by the coupling between spin-waves and the EM field 
in the FI. We have also studied the Fiske resonance in the S/I/FI/S junction and found that 
additional resonant structures appear due to the coupling between the ac Josephson current 
and the EM field in the I layer. 
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Appendix A: Derivation of Eq. (8) 

We integrate Eq. (1) with a narrow stripe A with infinitesimal width (dy) in the yz-plane, 

/ dS^-E z (y,t) = - f dS^-[fi H x (y,z,t) + M x (y,t)}, (A-l) 

J a dy J a 9t 

where E z (y,t) is confined to the ferromagnetic layer due to vanishing E z (y,t) in the S. Inte- 
grating Eq. (A-l) with respect to y and introducing the London penetration depth Al defined 
by 

l r ±0 ° 

Al = TT7 tt / dzH x (y,z,t), d F i « A L , (A-2) 

Eq. (A-l) becomes 

-d FI —E z (y,t) = fi ^H x (y,t)(d FI + 2\ L ) + d FI —M x (y,t). (A-3) 

In the same way, we integrate Eq. (2) over the cross-section area S' of the junction in the 
xz-plane and obtain 

H x (y,t) = ±-I M (y,t), I M (y,t)= I dSJ^{y,t), (A-4) 
w Js> 

where W is the width along the x-axis of the junction and Iyi(y,t) is the current at position 
y in the superconducting electrode. Substituting Eq. (A-4) into Eq. (A-3) and Eq. (4), we 
obtain the partial differential equation 

^V(y,t) = -^(d F1 + 2X h )p M (y,t) + ^ J™ dy'dt' X (y -y',t- t')I M (y> ,t>), (A5) 

where V(y,t) = d F E z (y,t) is the voltage across the FI. The Ampere's law in the FI is 

— I M (y,t) = -^-W-V(y,t) - WJ z 3 (y,t) - WJ^(y,t) (A6) 

Differentiating partially Eq. (A-5) with respect to y, Eq. (A-5) becomes 

where 

d f°° 

F(y,t) = ^J dy'dt'x(y -y',t- t')l M (y',t'). (A-8) 

To obtain the voltage equation coupled with spin waves in the FI, we transform the function 
F(y, t). When we execute the Fourier transformation on x(y ~y\t — t') and Iu.(y' t') referring 
to y 1 , F(y,t) is given by 

1 f°° 

F(y, t) = ^J dqdt'e^xiq, t - t')iqI M (q, t'). (A-9) 
Making use of the Fourier transformation of y in Eq. (A-6), 

iqlu(q, t) = -C'§-V(q, t) - WJ^q, t) - W J z Q (q, t), (A-10) 
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and substituting Eq. (A-10) into Eq. (A-9), F(y,t) is given by 



/oo 
dy'dt' X (y-y',t-t') 
-oo 



C-^V(y',t') + WJKy',t>) + WJ^y',t') 



(A-ll) 



Where C = epiW/dpi- Substituting Eq. (A-ll) into Eq. (A-7) and using the relation V(y,t) 
d F \E z {y,t), we have the partial differential equation 
d 2 V{y,t) 1 ra2i 



dy 2 c? 



fi L 



+ r FI gZ|^) + r FI dFI -1 r dy'dt' X (y -y',t- t')V(y',t' 
at dpi + 2Al mo j_oo 

+ — I— Jf(y,*) + — J dlL -- — r dy'dt'x{y-y',t-t r )J^(y',i(^-l2) 

Appendix B: Calculation of factor g n 

To obtain the complex coefficient g n , substituting Eq. (13) and Jj « J c sin(u;ji — fcny) into 
Eq. (12), we obtain the equation for g n , 

oo 2 

Que™ 3 * cos (k n y) [-u 2 n + //(fc n , -cuj)^ 2 - iT Flt i(k n , - Wj ) Wj ] = ^fM^H, -wjje^*"^) 

n=l A J 

where w n = (cpi7r/L)n and //(</, — wj) = 1 + x(g, — ajj)dpi/[(dpi + 2\j)fi ] and g denotes &h or 
k n . For both sides of Eq. (B-l), multiplying cos(k n y) and performing integration with respect 
to y from to L, we get g n as follows: 

9n = ~ ^/^(^H,-Wj)— B " 7T 7, (B-2) 

2 /" L 

B n = Y / dy cos (k n y) cos (k H y) , (B-3) 
^ JO 

2 r L 

Cfl = Lj d y cos ( k ™y) s ' m ( k ny) ■ ( B4 ) 

Appendix C: Linearized solution of LLG equation 

We outline the derivation of the magnetic susceptibility in the FI by using the LLG 
equation. We consider the situation, in which the direction of magnetization is perpendicular 
to the junction. In the experimental measurement of the Fiske resonance, the dc magnetic 
field is smaller than several tens of gauss. Therefore, we can neglect the gradient of the 
magnetization due to the applied dc magnetic field. Since the precessional angle of spin in 
the FI is usually very small even at the magnetic resonance, we linearize the LLG equation 
as follows: 

^ = - 7 M y (H K - M z j m ) + 7 M z Ad 2 M y - a^L, (Gl) 
^ = 7 M X (H K - M z /m) - -M : h, - jM z Ad 2 M x + a^, (C-2) 
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where Hk and M z /fiQ are an anisotropic field and a demagnetization field in the FI, respec- 
tively. We assume the solutions of Eq. (CT) and (C-2) of the form 

M^xifo.wj^e^'e-*™, i = x,y, (C-3) 

where Xi(g, coj) is a magnetic susceptibility of FI and h x is the amplitude of x-component of 
ac magnetic field and q is a wave number of spin wave. Substituting Eq. (C-3) into Eq. (C-l) 
and (C-2), the solutions of the linearized LLG equation are 

M x = xx(q,uj)h x , (C-4) 

Xx(g,wj) = 7M z -2 2 , (C-5) 

iZg — (1 + a JWj + zzaiis^j 

M y = Xy(q,uj)h x , (C-6) 

X s (?,wj) = -7^^ n T 2T2 T o o ' ( C ' 7 ) 

i2g — (1 + a z )ujj + izaiiswj 

ft s = ^B + V, n B = 1 (H K -M z / f i ), (C-8) 

where -£Tk is an anisotropic field in the FI and 77 is a stiffness of spin waves. 

Appendix D: Derivation of dc Josephson current in S/I/FI/S junction 

In this Appendix, we consider the Josephson junction including a nonmagnetic insulator 
(I) between the S and the FI, that is the S/I/FI/S junction. 

To calculate the Josephson current in the S/I/FI/S junction, we solve the Maxwell equa- 
tion as shown in the section II to add to equations describing the I as follows: 

rot [Ei(y,t)e z ]=-no§i [H l x (y,z,t)\ e x , (D-l) 

rot [Hl(y,t)e x ] - ej- [E l z {y,t)e z ] = J|(y, t)e g + <nE l z {y, t)e z , (D-2) 

where the subscript I in the above equations indicates fields in the non-magnetic insulator. e\ 
and o"i are a dielectric constant and a conductivity in the I, respectively. Jj is a Josephson 
current. In the same manner as the calculation in the section II, we integrate the Faraday's 
law with respect to the plane parallel to the yz-plane and the Ampere's law with respect to 
the plane parallel to the xz-plane. As a result, Ampere's and Faraday's laws are given by 
dEV-(y,t) , dEl(y,t) , dH^(y,t) . , XN dMJt) . 

+ fXo dHl(y,t) ((fa + AL)> (d . 3) 

H x (y,t) = ±li i = FI or I. (D-4) 

The definition of is same as that of Eq. (A-4). Performing the procedure of calculation to 
obtain Eq. (8) for Eqs. (D-3) and (D-4), we obtain equations describing ac voltages induced 
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by the ac Josephson current in the I and the FI layers as, 

d 2 V\y,t) 1 d 2 V l {y,t) l r dV l (y,t) h 1 dJ 3 (y,t) 

+ -9- 1 - 1 — h 



dy 2 



d 2 t 



dt 2e\ 2 J c 



dt 



di + \ L 



-eiAtO) 



A7 



d 2 V M (y,t) 
dy 2 



1 d 2 V FI {y,t) 



2e 1 
— fi (di + X L ), Fi = ——, 
n ei-Ki 



(D-5) 
(D-6) 



FI 



1 



+ — Tfi 



d 2 t 

dV Fl 



+ 



dm 



c|i d Fi + A L 9t y_ 



OC 



/oo 
-00 



at 



+ 



,-2 



2e Apj J c 



9Jj{y,t) h 



—00 
l 



9t 



2e Apj J c cZfi + Al dt 



dt ; X (t-t ; )Jj(y,t'), (D-7) 



<-FI — 3 CFI^O, "pi 

«FI 

where cZi is the thickness of I. ei and R\ are a dielectric constant and resistance per unit 
length in the I. Here, we adopt the iterative calculation to obtain the dc Josephson current 
due to the resonance between the electromagnetic field and Josephson current as shown in 
previous section. First, we calculate the voltage induced by the ac Josephson current in each 
layers. Josephson current is characterized by the phase difference between two S's. Therefore, 
we assume that the Josephson current flowing layers is given by Jj(y,t) = J c sin(wjt — kny), 
where kn = 2ir/j,oH ex (dm + di)/$o an d ojj is the Josephson frequency. For the boundary 
condition, we adopt the open-ended boundary condition for the reason mentioned in section 
II. In this case, the voltage expression satisfying the open-ended boundary condition is given 
by 



a; 



2e , 1 

— fi (dm + Al) , 1 fi = — ^— ■ 
h e F ii?Fl 



(D-8) 



V\y,t) = lm 



cos (k n y) 



,n=l 



j = I or FI, 



(D-9) 



where v 3 n are complex numbers and k n = nir/L. Substituting Eq. (D-9) into Eqs. (D-5) and 
(D-7), we can obtain the voltage in I and F layers as follows: 



V\y,t) 



h 

~e~L 



2 00 

E 

n=l 



-TuojBn + coj [uJ 2 - {uj l n ) 2 } C n 



H — — cos (ujt) 



{u\ - u 3 y + (riwj) 



h 2 -K) 2 ] + (W 

cos(k n y), 



sin (ojjt) 



(D-10) 



h ( cm 



v ^- t) = -7L{f 



2 00 

E 

n=l 



Mg(uj) sin(a;j) + cos(wjt))] cos{k n y), (D-ll) 



C \ 

Aj ) ^ WJ ^2 - fj,(-uj)uj + ir F iAt(-a;j) ' 



(D-12) 



where u\ and w n are given by (ci7r/L)n and (cfitt / L)n, respectively. B n and C n are same 
expressions as Eq. (15) and Eq. (16), respectively. 
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Next we calculate the dc Josephson current coupled with spin waves as a function of the 
dc voltage and of the external magnetic field. Since we consider y*( FI ) (y, t) as a perturbation, 
we can expand the sine function in terms of V l ( Fl \y, t). Within the first order term with 
respect to V 1( - Fl ^ (y,t), the dc Josephson current is approximately given by 

J dc « lim - / dt- dy,J c cos (ujt-kny)— / dt [V\y, t) + V F \y, t)] . (D-13) 

Substituting Eqs. (D-10) and (D-ll) into Eq. (D-13), we can obtain the analytic formula of 
the dc Josephson current in the S/I/FI/S junction as follows: 

J*, = & + J™, (D14) 



4^£ „,„ Xrr (D ' 15) 

4A I n=l L( w n) - W jJ + ( r i^j) 

4? = ^£*^), (D-16) 

4A FI n=l 

_ [/A^J^fi^j + fi"(ujj)Loj] + fi"(u}j) [u>n — [i'(ujj)ujj + ^"(ajj)rFiu;j] 

[a;2 - //'(wj)^ + n"(uj)r F iuj] 2 + [^'(u;j)rFiu;j + 2 

(D-17) 

where and w n are given by (ci7r/L)n and (cfjtt / L)n, respectively. F n ((f>) is same function as 
Eq. (25). Jjj° comes from the resonance between the electromagnetic field generated inside the 
I and the ac Josephson current. Jjpj originates in the resonance between the electromagnetic 
field generated inside the FI and the ac Josephson current. 



15/16 



J. Phys. Soc. Jpn. Full Paper 

References 

1) B.D. Josephson, Phys. Lett. 1 (1962) 251. 

2) M. D. Fiske, Rev. Mod. Phys. 36 (1964) 221. 

3) R. E. Eck, D. J. Scalapino, and B. N. Taylor, Phys. Rev. Lett. 13 (1964) 15. 

4) D. D. Coon and M.D Fiske, Phys. Rev. 138 (1965) A744. 

5) I. O. Kulik, JETP Lett. 2 (1965) 84. 

6) A. Barone and G. Paterno, Physics and Applications of the Josephson Effect (Wiley, New York, 
1982). 

7) A. A. Golubov, M. Yu. Kupriyanov, and E. Ilfichev, Rev. Mod. Phys. 76 (2004) 411. 

8) A. I. Buzdin, Rev. Mod. Phys. 77 (2005) 935. 

9) V. V. Ryazanov, V. A. Oboznov, A. Yu. Rusanov, A. V. Vcrctcnnikov, A. A. Golubov, and J. Aarts, 
Phys. Rev. Lett. 86 (2001) 2427. 

10) T. Kontos, M. Aprili, J. Lesueur, F. Genet, B. Stcphanidis, and R. Boursier, Phys. Rev. Lett. 89 
(2002) 137007. 

11) Z. Nussinov, A. Shnirman, D. P. Arovas, A. V. Balatsky, and J. X. Zhu, Phys. Rev. B 71 (2005) 
214520. 

12) S. Takahashi, S. Hikino, M. Mori, J. Martinck and S. Maekawa, Phys. Rev. Lett. 99 (2007) 057003. 

13) C. Bell, S. Milikisyants, M. Hubcr, and J. Aarts, Phys. Rev. Lett. 100 (2008) 047002. 

14) M. Houzet, Phys. Rev. Lett. 101 (2008) 057009. 

15) S. Hikino, M. Mori, S. Takahashi, and S. Maekawa, J. Phys. Soc. Jpn. 77 (2008) 053707. 

16) F. Konschelle and A. Buzdin, Phys. Rev. Lett. 102 (2009) 017001. 

17) T. Yokoyama and Y. Tscrkovnyak, Phys. Rev. B 80 (2009) 104416. 

18) I. Pctkovic, M. Aprili, S. E. Barnes, F. Beuneu, and S. Maekawa, Phys. Rev. B 80 (2009) 220502. 

19) A. F. Volkov and K. B. Efctov, Phys. Rev. Lett. 103 (2009) 037003. 

20) M. Aprili, Spin and Phase dynamics in Ferromagnetic Josephson junctions, KN-3-2, "9 th 
International Conference on Materials and Mechanisms of Superconductivity^ , Tokyo Japan. 

21) S. Kawabata, Y. Asano, Y. Tanaka, A. A. Golubov, and S. Kashiwaya, Phys. Rev. Lett. 104 (2010) 
117002. 

22) S. Kawabata, S. Kashiwaya, Y. Asano, Y. Tanaka, and A. Golubov, Phys. Rev. B 74 (2006) 180502 
(R). 

23) B. Hillebrands and K. Ounadjela, Spin dynamics in confined magnetic structures II (Springer- 
Verlag Berlin Heidelberg, New York, 2003) 

24) Y. Kajiwara, K. Harii, S. Takahashi, J. Ohe, K. Uchida, M. Mizuguchi, H. Umezawa, H. Kawai, 
K. Ando, K. Takanashi, S. Maekawa, and E. Saitoh Nature 464 (2010) 262. 

25) D. D. Stancil and A. Prabhakar, Spin Waves Theory and Applications (Springer Science + Business 
Media, LLC 2009) 

26) M. P. Lisitskiy and M. V. Fistul, Phys. Rev. B 81 (2010) 184505. 

27) For L = 1 mm and cpi = 10 7 m/s, the magnitude of dc voltage (Vb) is about tens of fiV. 

28) M. Pajda, J. Kudrnovsky, I. Turek, V. Drchal, and P. Bruno, Phys. Rev. B 64 (2001) 174402. 

29) U. Ebels, P. E. Wigen, and K. Ounadjela, Europhys. Lett. 46 (1999) 94. 



16/16 



